Electron- and neutrino-nucleus scattering in the impulse approximation regime 
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A quantitative understanding of the weak nuclear response is a prerequisite for the analyses of 
neutrino experiments such as K2K and MiniBOONE, which measure energy and angle of the muons 
produced in neutrino-nucleus interactions in the energy range 0.5 -- 3 GeV and reconstruct the 
incident neutrino energy to determine neutrino oscillations. In this paper we discuss theoretical 
calculations of electron- and neutrino-nucleus scattering, carried out within the impulse approxi- 
mation scheme using realistic nuclear spectral functions. Comparison between electron scattering 
data and the calculated inclusive cross section off oxygen, at beam energies ranging between 700 
and 1200 MeV, show that the Fermi gas model, widely used in the analysis of neutrino oscillation 
experiments, fails to provide a satisfactory description of the measured cross sections, and inclusion 
of nuclear dynamics is needed. 
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I. INTRODUCTION 

The field of neutrino physics is rapidly developing af- 
ter atmospheric neutrino oscillations and solar neutrino 
oscillations have been established 0,0,11. Recently, the 
SK Collaboration has found evidence of the oscillatory 
signature in atmospheric neutrinos, improving the deter- 
mination of Am^ U , and K2K experiment has confirmed 
the oscillations of atmospheric neutrinos at 99.995% CL 
1^ 1^ . These neutrino experiments measure energy and 
angle of muons produced in neutrino-nucleus interactions 
and reconstruct the incident neutrino energy, which de- 
termines the neutrino oscillations. K2K took data in the 
El, = 0.5 — 3 GeV region, and the recent L/E analy- 
sis of the SK atmospheric neutrinos is mainly based on 
the dataset extending from 0.5 to 25 GeV. JPARC and 
NuMI neutrino experiments 0, ll| propose to measure 
Ufj_ — > Ue oscillations and determine Am^ with 1% ac- 
curacy and sin^ 20i3 above 0.006, using a narrow-band 
neutrino beam at E,, = 0.8 GeV (JPARC) and 2.0 GeV 
(NuMI, off-axis). 

In view of these developments, it is vital that theoret- 
ical calculations of cross sections and spectra achieve an 
accuracy comparable to the experimental one, which in 
turn requires that the nuclear response to weak interac- 
tions be under control at a quantitative level. 

At E^=3 GeV or less, quasi-elastic scattering and 
quasi-free A production are the dominant neutrino- 
nucleus processes. However, reactions in this energy 



regime are associated with a wide range of momentum 
transfer, thus involving different aspects of nuclear struc- 
ture. 

Four decades of electron-nucleus scattering experi- 
ments have unequivocally shown that the mean-field ap- 
proximation, underlying the nuclear shell model, does 
not provide a fully quantitative account of the data (see, 
e.g., Ref. |2| and references therein). When the mo- 
mentum transfer involved is large, dynamical nucleon- 
nucleon (NN) correlations are known to be important, 
and a description of nuclear structure beyond the mean- 
field picture is needed. On the other hand, neutrino- 
nucleus reactions also occur, in fact rather appreciably, 
with a small momentum transfer. Comparison between 
the data at < 0.2 GeV'^ and the predictions of the 
Fermi gas (FG) model jT3|, showing a sizable deficit of 
events [y, [ll| , suggests that a more realistic description 
of both nuclear properties and the reaction mechanism is 
indeed required. 

In this paper we discuss the extension of the many- 
body theory of electron- nucleus scattering (see, e.g., Ref. 
[l^ and references therein) to the case of neutrino- 
induced reactions. We focus on the energy range 0.7—1.2 
GeV and analyze inclusive scattering of both electrons 
and neutrinos off oxygen, the main target nucleus in SK, 
K2K and other experiments. The quasi-elastic and quasi- 
free A production cross sections obtained from the FG 
model [io, [l3 are compared to the results of the many- 
body approach developed in Ref. Q, extensively used 
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to analyze electron scattering data at beam energy up to 
few GeV El El El. Preliminary versions of the mate- 
rials in this paper have appeared in the Proceedings of 
Nulnt04 '^,17', 'T^. 

Section ^ is devoted to a summary of the formalism 
employed to calculate the electron-nucleus cross section 
at high momentum transfer, as well as to the discussion 
of the main ingredients entering its definition: the nu- 
clear spectral function, the elementary cross section de- 
scribing electron scattering off a bound nucleon and the 
folding function embodying the main effects of final state 
interactions. 

In Section imi the results of our approach are compared 
to inclusive electron scattering data at 0.2 ^ ^ 0.6 
GeV^, while in Section Hvl we outline the extension of 
the formalism to the case of charged current neutrino- 
nucleus scettering. Finally, our conclusions are stated in 
Section V. 



II. MANY-BODY THEORY OF THE 
ELECTROWEAK NUCLEAR RESPONSE 

A. Electron-nucleus cross section 



The differential cross section of the process 



A 



X 



(1) 



in which an electron carrying initial four-momentum 
k = [Ef. , k) scatters off a nuclear target to a state of 
four-momentum k' = (£'g,k'), the target final state be- 
ing undetected, can be written in Born approximation as 
(see, e.g., Ref. ^]) 

where a is the fine structure constant and = —q^ = 
— j/^, q = k — k' = (t^, q) being the four momentum 
transfer. 

The leptonic tensor, that can be written, neglecting 
the lepton mass, as 

(3) 



2 yt^k^ ~\- k^k^ {kk')] 



is completely determined by electron kinematics, whereas 
the nuclear tensor W^'^ contains all the information on 
target structure. Its definition involves the initial and 
final hadronic states |0) and \X), carrying four-momenta 
po and pxj respectively, as well as the nuclear electro- 
magnetic current operator J'^: 

W^'' ^Y.'^Q\J^^\^){X\.rm^^\pa + q-px) , (4) 

X 

where the sum includes all hadronic final states. 

Calculations of W^'^ at moderate momentum transfers 
(|q| < 0.5 GeV) can be carried out within nuclear many- 
body theory (NMBT), using nonrelativistic wave func- 
tions to describe the initial and final states and expand- 
ing the current operator in powers of jql/m, m being the 



nucleon mass (see, e.g., Ref. Hy). On the other hand, 
at higher values of |q|, corresponding to beam energies 
larger than ^ 1 GeV, the description of the final states 
\X) in terms of nonrelativistic nucleons is no longer pos- 
sible. Galculations of W^" in this regime require a set 
of simplifying assumptions, allowing one to take into ac- 
count the relativistic motion of final state particles car- 
rying momenta ~ q as well as the occurrence of inelastic 
processes, leading to the appearance of hadrons other 
than protons and neutrons. 



B. The impulse approximation 

The main assumptions underlying the impulse approx- 
imation (lA) scheme are that i) as the spatial resolu- 
tion of a probe delivering momentum q is ~ l/|q|, at 
large enough |q| the target nucleus is seen by the probe 
as a collection of individual nucleons and ii) the parti- 
cles produced at the interaction vertex and the recoiling 
(A — I)-nucleon system evolve indipendently of one an- 
other, which amounts to neglecting both statistical corre- 
lations due to Pauli blocking and dynamical Final State 
Interactions (FSI), i.e. rescattering processes driven by 
strong interactions. 

In the lA regime the scattering process off a nuclear 
target reduces to the incoherent sum of elementary pro- 
cesses involving only one nucleon, as schematically illus- 
trated in Fig. nj 




FIG. 1: (Color online) Pictorial representation of the lA 
scheme, in which the nuclear cross section is replaced by the 
incoherent sum of cross sections describing scattering off indi- 
vidual bound nucleons, the recoiling (A — l)-nucleon system 
acting as a spectator. 



Within this picture, the nuclear current can be written 
as a sum of one-body currents 



(5) 



while the final state reduces to the direct product of the 
hadronic state produced at the electromagnetic vertex, 
carrying momentum p^. and the {A — I)-nucleon residual 
system, carrying momentum p-ji = q — (for simplicity, 
we omit spin indices) 



\x) \x,p^) (g> |7^,PK) 



(6) 
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Using Eq. (jBJ we can rewrite the sum in Eq. replacing 
X ^dV|7^,PK)(P7^,7^| . (7) 



Substitution of Eqs. lO-l(7|) into Eq. Q and insertion of 
a complete set of free nucleon states, satisfying 

J d'p\N,p){p,N\=I , (8) 

results in the factorization of the current matrix clement 



IJ'^IX) 



:(0|7^,p7^;N,-p7^) 



y/p'^ + m2 
X ^(-pK,iV|jf|x,p,) , 



(9) 



leading to 

W^'' = f dVd'Px|(0|7^,P7^;N,-p7^)p 



X ^ (-pK,N|jf|a:,p,)(p,,x|jT|N,-pK) (10) 

i 

X S^^\q-pn-Px)S{u + Eo-En-E,), 
where Ep^ = \/\ptz\^ + m? ■ Finally, using the identity 

5{v + Eo-En-E,) = J dE6{E -m + Eo-En) 
X 5{v~E + m~E^) , (11) 
and defining the target spectral function as 0| 

P{p,E) = ^|(0|7^,-p;N,p)|2 
n 

xS{E-m + Eo-En) , (12) 
we can rewrite Eq. I^J in the form 



E^ 



P{p,E), (13) 



with Ep = \/\p'^\ + m? and 

= E(P,N|jrk,P + q)(p + q,x|jnN,p) 



X 5{l' + ^p2 + TO^ — i?a; 



(14) 



The quantity defined in the above equation is the ten- 
sor describing electromagnetic interactions of the i-th nu- 
cleon in free space. Hence, Eq. (|14|l shows that in the lA 



scheme the effect of nuclear binding of the struck nucleon 
is accounted for by the replacement 



with (see Eqs. ^ and lO) 

V = Ex — \/p'^ +m? 



(15) 



= v + E^-En- Vp^T^ 
= v — E + m — \/p^ + m^ , 



(16) 



in the argument of iff. This procedure essentially 
amounts to assuming that: i) a fraction 5v of the en- 
ergy transfer goes into excitation energy of the spectator 
system and ii) the elementary scattering process can be 
described as if it took place in free space with energy 
transfer v — v — 5v. This interpretation emerges most 
naturally in the |p| <C m limit, in which Eq. (|15|l yields 
5v = E. 

Collecting together all the above results we can finally 
rewrite the doubly differential nuclear cross section in the 
form 



daiA 

dn 



+ {A~Z) 



^ dcTep 



dap. 



dilp'dEpi 



dfle'dEg' 
5{v-E + m~E.,), (17) 



where dcFej^ /dfle'dEe' {N = n,p denotes a proton or a 
neutron) is the cross section describing the elementary 
scattering process 



3(A:) + N{p) 



e'{k') + x{p + q) 



given by 



da, 



eN 



o? E' m 



drie'dE,, Q*eI Ep^'"''"^'' ' 



(18) 



(19) 



stripped of both the flux factor and the energy conserving 
(5-function. 



C. The nuclear spectral function 

In NMBT the nucleus is seen as a system of A nucle- 
ons whose dynamics are described by the nonrelativistic 
hamiltonian 



A 2 -4 
i=l j>i=l k>j>i=l 



E ' 



(20) 



where pi is the momentum of the i-th nucleon, while Vij 
and Vijk are two- and three-nuclcon interaction poten- 
tials, respectively. 

The two-nucleon potential, that reduces to the Yukawa 
one-pion-exchange potential at large internucleon dis- 
tance, is obtained from an accurate fit to the available 



4 



data on the two-nucleon system, i.e. deuteron proper- 
ties and ~ 4000 NN scattering phase shifts at energies 
up to the pion production threshold j^^] . The additional 
three-body term Vijk has to be included in order to ac- 
count for the binding energies of the three-nucleon bound 
states j23i] . 

The many-body Schrodinger equation associated with 
the Hamiltonian of Eq. H20|l can be solved exactly, using 
stochastic methods, for nuclei with mass number A < 10. 
The energies of the ground and low-lying excited states 
are in excellent agreement with the experimental data 
[2^. Accurate calculations can also be carried out for 
uniform nucleon matter, exploiting translational invari- 
ace and using either a variational approach based on clus- 
ter expansion and chain summation techniques |2^, or 
G-matrix perturbation theory (2^. 

Nonrelativistic NMBT provides a fully consistent com- 
putational framework that has been employed to obtain 
the spectral functions of the few- nucleon systems, ]ia,y- 
ing A= 3 0, m m and 4 [H US IP, as well as of 
nuclear matter, i.e. in the limit A ^ 00 with Z=A/2 
[3^ Ib^ I . Calculations based on G-matrix perturbation 
theory have also been carried out for oxygen [s^ . Issf . 

The spectral functions of different nuclei, ranging from 
Carbon to Gold, have been modeled using the Local Den- 
sity Approximation (LDA) 'l6l , in which the experimen- 
tal information obtained from nucleon knock-out mea- 
surements is combined with the results of theoretical cal- 
culations of the nuclear matter P(p, E) at different den- 
sities 16]. 

Nucleon removal from shell model states has been ex- 
tensively studied by coincidence (e, e'p) experiments (see, 
e.g., Ref. 01). The corresponding measured spectral 
function is usually written in the factorized form 



Pmf{p, E)=Y1 ^« \MP)\^Fn{E - En) 



(21) 



where (^„(p) is the momentum-space wave function of the 
single particle shell mode state n, whose energy width is 
described by the function Fn{E — En) ■ The normalization 
of the n-th state is given by the so called spectroscopic 
factor Zn < 1, and the sum in Eq. (|21|l is extended to 
all occupied states of the Fermi sea. Hence, Pmf{p,E) 
vanishes at |p| larger than the Fermi momentum pf ^ 
250 MeV. Note that in absence of NN correlations Fn{E- 
En) shrinks to a (5-function, Z„ = 1 and Eq. (|21|l can be 
identified with the full spectral function. 

Strong dynamical NN correlations give rise to virtual 
scattering processes leading to the excitation of the par- 
ticipating nucleons to states of energy larger than the 
Fermi energy, thus depleting the shell model states within 
the Fermi sea. As a consequence, the spectral function 
associated with nucleons belonging to correlated pairs ex- 
tends to the region of |p| ^ pf and E ^ cf, where cf 
denotes the Fermi energy, typically < 30 MeV. 

The correlation contribution to P(p, E) of uniform nu- 
clear matter has been calculated by Benhar et al for 
a wide range of density values Within the LDA 



scheme, the results of Ref. can be used to obtain 
the corresponding quantity for a finite nucleus of mass 
number A from 

P,orr{p,E)^ J d\ pa{v)P^J;.{p,E-p^ pa{v)) , (22) 

where pa{i^) is the nuclear density distribution and 
P(^rr {Pt E; p) is the correlation component of the spec- 
tral function of uniform nuclear matter at density p. 

Finally, the full LDA nuclear spectral function can be 
written 



PlDa{P, E) = Pmf{P, E) + PcorriP. E) , 



(23) 



the spectroscopic factors Z„ of Eq. (|21|l being con- 
strained by the normalization requirement 



j d^p dE Plda{p,E) = 1 



(24) 
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FIG. 2: (Color online) Three-dimensional plot (left panel) 
and scatter plot (right panel) of the oxygen spectral function 
obtained using the LDA approximation described in the text. 

The LDA spectral function of ^^O obtained combin- 
ing the nuclear matter results of Ref. and the Saclay 
(e, e'p) data [S^l is shown in Fig. |21 The shell model con- 
tribution Pmf{p, E) accounts for 80 % of its normal- 
ization, whereas the remaining ^ 20 % of the strength, 
accounted for by Pcorrip, E), is located at high momen- 
tum (IpI 3> pf) and large removal energy (E ^ e^)- 
It has to be emphasized that large E and large p are 
strongly correlated. For example, ~ 50 % of the strength 
at IpI = 320 MeV is located at > 80 MeV. 

The LDA scheme rests on the premise that short range 
nuclear dynamics is unaffected by surface and shell ef- 
fects. The validity of this assumption is confirmed by 
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theoretical calculations of the nucleon momentum distri- 
bution, defined as 

n{p) = J dE P{p,E) (25) 
= (Ola^flplO) , (26) 

where and Op denote the creation and annihilation op- 
erators of a nucleon of momentum p. The results clearly 
show that for A> 4 the quantity n{p)/A becomes nearly 
independent of A in the region of large |p| (> 300 MeV), 
where NN correlations dominate (see, e.g., Ref. [Ill)- 

In Fig. 13 the nucleon momentum distribution of ^^O, 
obtained from Eq. H25|l using the LDA spectral function 
of Fig. 121 is compared to the one resulting from a Monte 
Carlo calculation ^3 , carried out using the definition of 
Eg. H26|l and a highly realistic many-body wave function 
|4l|. For reference, the FG model momemtum distribu- 
tion corresponding to Fermi momentum pp ~ 221 MeV 
is also shown by the dashed line. It clearly appears that 
the fT.(p) obtained from the spectral function is close to 
that of Ref. ^3 > while the FG distribution exhibits a com- 
pletely different behaviour. 




12 3 4 



P [fm-'] 

FIG. 3: (Color online) Momentum distribution of nucleons 
in the oxygen ground state. Solid line: LDA approximation. 
Dashed line: FG model with Fermi momentum pp — 221 
MeV. Diamonds: Monte Carlo calculation carried out by S.C. 
Pieper ^) using the wave function of Ref. |4l| . 



A direct measuremet of the correlation component of 
the spectral function of ^^C, obtained measuring the 
(e, e'p) cross section at missing momentum and energy 
up to ~ 800 MeV and ^ 200 MeV, respectively, has been 
recently carried out at Jefferson Lab by the E97-006 Col- 
laboration .3^ . The data resulting from the preliminary 
analysis appear to be consistent with the theoretical pre- 
dictions based on LDA. 



D. Final state interactions 

The occurrence of strong FSI in electron nucleus scat- 
tering has long been experimentally established. The re- 
sults of a number of (e, e'p) measurements covering the 
kinematical domain corresponding to 0.5 iS 1z 8.0 
GeV2 l42l|4il3E3, clearly show that the flux of out- 
going protons is strongly suppressed, with respect to the 
lA predictions. The observed attenuation ranges from 
20-40 % in Carbon to 50-70 % in Gold. 

The inclusive (e, e') cross section, being only sensitive 
to rescattering processes taking place within a distance 
~ l/|q| of the electromagnetic vertex, is obviously much 
less sensitive to FSI than the coincidence {e,e'p) cross 
section. The latter is in fact affected by rescatterings 
occurring over the distance ^ Ra, Ra being the nuclear 
radius, travelled by the struck particle on its way out 
of the target. However, FSI effects become appreciable, 
indeed dominant, in the low v region, where the inclusive 
cross section is most sensitive to the high momentum and 
high removal energy tails of the nuclear spectral function. 

In quasi-elastic inclusive processes FSI produce two ef- 
fects: i) an energy shift of the cross section, due to the 
fact that the struck nucleon feels the mean field gener- 
ated by the spectator particles and ii) a redistribution 
of the strength, leading to the quenching of the peak 
and the enhancement of the tails, to be ascribed to the 
occurrence of NN rescattering processes that couple the 
one particle-one hole final state to more complicated n 
particle-n holes configurations. 

Early attempts to include FSI effects were based on the 
optical potential model j4^ . However, while providing a 
computationally practical scheme to account for the loss 
of flux in the one-nucleon removal channel, this model 
relies on the mean field picture of the nucleus, and does 
not include the effect of dynamical NN correlations. 

A different approach, based on NMBT and a gener- 
alization of Glauber theory of high energy proton scat- 
tering ^3 has been proposed by Benhar et al. |0| in 
the early 90s. This treatment of FSI, generally referred 
to as Correlated Glauber Approximation (CGA) rests on 
the assumptions that i) the struck nucleon moves along 
a straight trajectory with constant velocity (eikonal ap- 
proximation), and ii) the spectator nucleons are seen by 
the struck particle as a collection of fixed scattering cen- 
ters (frozen approximation). 

Under the above assumptions the expectation value of 
the propagator of the struck nucleon in the target ground 
state can be written in the factorizcd form 

f/p+q(i) = K+^m^^lit) , (27) 

where C^p+q is the free space propagator, while FSI effects 
are described by the quantity (i? = (ri, . . . , r^i) specifies 
the target configuration) 

u^^iit)^m^^i{R;tm , (28) 
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with 



ui 



,(|rij+vt'|) 



(29) 



In Eq. 129|) Vij = r,; —Vj and Wp+q(|r|) is the coordinate- 
space NN scattering t- matrix at incident momentum p + 
q, usually parametrized in terms of total cross section, 
slope and real to imaginary part ratio. At large |q|, p + 
q « q and the eikonal propagator of Eq. (|28|l becomes a 
function of t and the momentum transfer only. 

Note that U^^^{R]t) is simply related to the nuclear 
transparency Ta , measured in coincidence (e, e'p) exper- 
iments millEllii, through 



Ta = lim(0||[/^^^(i?;i)|^|0) . 



(30) 



The results displayed in Fig. 0] |[4^] show that both the 
magnitude and the A- and Q^-dependence of the trans- 
parencies of Carbon, Iron and Gold obtained from the 
approach of Ref. [iJI and LDA are in good agreement 
with the experimental data. Note that at low FSI 
lead to a 20 (40) % effect in Carbon (Iron). Neglect- 
ing this effect, i.e. setting Ta{Q'^) = 1, would be utterly 
incompatible with the data. 
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FIG. 4: (Color online) Q^-dependence of the transparency of 
Carbon, Iron and Gold ji^j . calculated using LDA and the 
approach of ref.|l4j|. The data points are taken from Refs.|42j| 
(crosses), (diamonds) and 0,^3 (squares). Note that, 
in absence of FSI, Ta(Q^) = 1. 

From Eqs. (|^ and it follows that within the ap- 
proach of Ref. 14] the energy shift and the redistribution 
of the inclusive strength are driven by the real and the 
imaginary part of the NN scattering amplitude, respec- 
tively. At large q the imaginary part of Wq, correspond- 
ing to the real part of U^^^ , is dominant. Neglecting 
the contribution of the real part of Wq altogether, the 
CCA quasi-elastic inclusive cross section can be written 



as a convolution integral, involving the cross section eval- 
uated within the lA, i.e. in absence of FSI, and a folding 
function embodying FSI effects: 



da 



dQ,pidv 



the folding function /q(t^) being defined as 
dt 



27r 



(32) 



The above equations clearly show that the strength of 
FSI is measured by both Ta and the width of the folding 
function. In absence of FSI Uq^^ {R;t) = 1, implying in 
turn T/i = 1 and fq{v) — S{h'). 

Dynamical NN correlations stronlgy affect the shape 
of the folding function of Eq. Due to the strong 

repulsive core of the NN force, the joint probability of 
finding two nucleons at positions and rj, driving the 
occurrence of rescattering processes in the final state, is 
strongly suppressed at jr^ — r^l < 1 fm. As a conse- 
quence, inclusion of correlation effects within the frame- 
work of NMBT leads to a strong quenching of FSI effects, 
with respect to the predictions of the independent parti- 
cle model. 

In principle, the real part of the NN scattering ampli- 
utde can be explicitely included in Eq. (|^ and treated 
on the same footing as the imaginary part. However, its 
effect turns out to be appreaciable only at i ^ 0, when 
the attenuation produced by the imaginary part is weak. 
The results of numerical calculations show that an ap- 
proximate treatment based on the use of a time indepen- 
dent optical potential is indeed adequate to describe the 
energy shift produced by the real part of Wq |16| , whose 
size of ~ 10 MeV is to be compared to a typical electron 
energy loss of few hundreds MeV. 



III. COMPARISON TO ELECTRON 
SCATTERING DATA 

We have employed the formalism described in the pre- 
viuos Sections to compute the inclusive electron scatter- 
ing cross section off oxygen at 0.2 ^ ^ 0.6 GeV^. 

The lA cross section has been obtained using the LDA 
spectral function shown in Fig. |21and the nucleon tensor 
defined by Eq. H14|l . that can be written as 



..N 



p 



, (33) 



where p = {Ep , k) and the off-shell four momentum 
transfer q is defined by Eqs. H15|) and (|16|l . The 
two structure functions and W2 are extracted from 
electron-proton and electron-deuteron scattering data. 



7 



In the case of quasi-elastic scattering they are simply re- 
lated to the electric and magnetic nucleon form factors, 
Gem and Gmm i through 



Am? 



1 



1 — /Am? 



3^ \ r2 

2m ' 



(34) 



2m 



Gl 



~2 

^ G^ 



Am^ 



(35) 



Numerical calculations have been carried out using the 
Hohler-Brash parameterization of the form factors [iol 
, resulting from a fit which includes the recent Jeffer- 
son Lab data 51]. 

In the kinematical region under discussion, inelastic 
processes, mainly quasi-free A resonance production, are 
also known to play a role. To include these contribu- 
tions in the calculation of the inclusive cross section, we 
have adopted the Bodek and Ritchie parametrization of 
the proton and neutron structure functions j52| | , covering 
both the resonance and deep inelastic region. 

The folding functions describing the effect of NN 
rescattering in the final state have been computed from 
Eq. (|^ with the eikonal propagator U^^^{R;t) ob- 
tained using the parametrization of the NN scattering 
amplitude of Ref. js^ and the medium modified NN 
cross sections of Ref. [s^. The integrations involved 
in Eq. (|28|l have been carried out using Monte Carlo 
configurations sampled from the probability distribution 
associated with the oxygen ground state wave function 
of Ref. 

The effect of the real part of the NN scattering am- 
plitude has been approximated including in the energy 
conserving ^-function of Eq. H17|l the real part of the 
optical potential felt by a nucleon of momentum p + q 
embedded in uniform nuclear matter at equilibrium den- 
sity. 

In Figs. I5I8I the results of our calculations are com- 
pared to the data of Ref. [s^, corresponding to beam 
energies 700, 880, 1080 and 1200 MeV and electron scat- 
tering angle 32°. For reference, the results of the FG 
model corresponding to Fermi momentum pi;- = 225 MeV 
and average removal energy e = 25 MeV are also shown. 

Overall, the approach described in the previous Sec- 
tions, involving no adjustable parameters, provides a 
fairly accurate account of the measured cross sections 
in the region of the quasi-free peak. On the other hand, 
the FG model, while yielding a reasonable description at 
beam energies 1080 and 1200 MeV, largely overestimates 
the data at lower energies. The discrepancy at the top 
of the quasi-elastic peak turns out to be ~ 25 % and ~ 
50 % at 880 and 700 MeV, respectively. 

The results of NMBT and FG model also turn out to be 
sizably different in the dip region, on the right hand side 
of the quasi-elastic peak, while the discrepancies become 



'0(e,e') 

E = 700 MeV, 6=32° 
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FIG. 5: (Color online) Cross section of the process ^^0(e,e') 
at beam energy 700 MeV and electron scattering angle 32°. 
Solid line: full calculation, carried out within the approach 
described in SectionlH] Dot-dash line: lA calculation, carried 
out neglecting FSI effects. Dashed line: FG model with pf = 
225 MeV and e — 25 MeV. The experimental data are from 
Ref.iS^. 
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FIG. 6: (Color online) Same as in Fig. |^ but for beam energy 
880 MeV. 



less pronounced at the A-production peak. However, it 
clearly appears that, independent of the employed ap- 
proach and beam energy, theoretical results significantly 
underestimate the data at energy transfer larger than the 
pion production threshold. 

In view of the fact that the quasi-elastic peak is cor- 
rectly reproduced (within an accuracy of ~ 10 %), the 
failure of NMBT to reproduce the data at larger w may 
be ascribed to deficiencies in the description of the ele- 
mentary electron- nucleon cross section. In fact, as ilus- 
trated in Fig. 13 the calculation of the lA cross section 
at the quasi-elastic and A production peak involves inte- 
grations of P{p, E) extending over regions of the (p, E) 
plane almost exaclty overlapping one another. 

To gauge the uncertainty associated with the descrip- 
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FIG. 7: (Color online) Same as in Fig. |^ but for beam energy 
1080 MeV. 
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FIG. 8: (Color online) Same as in Fig. |^ but for beam energy 
1200 MeV. 



tion of the nucleon structure functions Wi and , 
we have compared the electron- proton cross sections ob- 
tained from the model of Ref. 153] to the ones obtained 
from the H2 model of Ref. [56] and from a global fit 
including recent Jefferson Lab data |^3|- The results of 
Fig. Unishow that at Ee = 1200 MeV and 9 = 32° the 
discrepancy between the different models is not large, 
being ^ 15 % at the A production peak. It has to be 
noticed, however, that the models of Refs. EE 
have all been obtained fitting data taken at electron beam 
energies larger than 2 GeV, so that their use in the kine- 
matical regime discussed in this work involves a degree 
of extrapolation. 

On the other hand, the results obtained using the ap- 
proach described in this paper and the nucleon structure 
functions of Ref. |52 are in excellent agreement with the 
measured (e, e') cross sections at beam energies of few 
GeV l()j. As an example, in Fig. E]we show a compar- 
ison between the calculated ^^C(e, e') cross section and 
the Jefferson Lab data of Ref. jSJl at i?e = 4 GeV and 
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FIG. 9: (Color online) The solid and dashed lines enclose 
the integration regions in the (p, E) plane relevant to the 
calculation of the lA cross section at the top of the quasi- 
elastic and A production peak, respectively, for beam energy 
1200 MeV and scattering angle 32°. 
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FIG. 10: (Color online) Cross section of the process e-^N ^ 
e'-f X above pion production threshold, at beam ener gy 1 200 
MeV and scattering angle 32° . Solid line: fit of Ref. for 
ep scattering; dashed line: fit of Ref. [H^ for ep scattering; 
diamonds: fit of Ref. for ep scattering; dot-dash line: fit 
of Ref. [s^l for en scattering. 



d = 30°. The corresponding FG result is also displayed, 
for reference. 

Fig. ^1 also shows the prediction of the Bodek and 
Ritchie fit for the neutron cross section, which turns out 
to be much smaller than the proton one. It is on account 
of this difference that we have chosen to adopt the fit of 
Ref. I as it allows for a consistent inclusion of proton 
and neutron contributions, both resonant and nonreson- 
sant, to the nuclear cross section. In this regard, it has 
to be pointed out that the nonresonant backround is not 
negligible. As illustrated in Fig. ^1 for beam energy 
1200 MeV and scattering angle 32° it provides ^ 25 % of 
the cross section at energy transfer corresponding to the 
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FIG. 11: (Color online) Cross section of the process e+^^C — > 
e' + X at beam energy 4 GeV and scattering angle 32° . The 
solid line has been obtained using the approach described in 
this work, while the dashed line shows the results of the FG 
model corresponding to pp = 221 MeV and e = 25 MeV. The 
experimental data are from Ref. [s^. 
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FIG. 12: (Color online) lA cross section of the process 
^^0{e,e') at beam energy 1200 MeV and scattering angle 
32°. Dashed line: quasi-elastic; dots: quasi-free A produc- 
tion; dashes: nonresonant background; solid line: total. The 
experimental data are from Ref. [5R| . 



A peak. 

The folding function described in Section ^ accounts 
for the FSI between a nucleon and the spectator system. 
For this reason, the results shown in Figs. I5I8I have been 
obtained folding with fq{v) only the quasielastic compo- 
nent of the lA cross section. Particles other than protons 
and neutrons, that can be produced at the electromag- 
netic vertex, also have FSI, but they are more difficult 
to describe. However, the inelastic part of the lA cross 
section, being rather smooth, is unlikely to be strongly 
affected by FSI. To gauge the possible relevance of ne- 
glecting FSI in the inelastic channels we have computed 
the cross section at incident energy 880 MeV and scat- 



tering angle 32° folding the total lA result. Comparison 
with the results displayed in Fig. shows a difference of 
^ 0.5 % at the top of the A-production peak. 



IV. CHARGED CURRENT 
NEUTRINO-NUCLEUS CROSS SECTION 

The Born approximation cross section of the weak 
charged current process 

iyi + A~^£- +X , (36) 

can be written in the form (compare to Eq. Q) 

da Ik' I 



dflidEi 327r2 Ikl 



(37) 



where G = GpcosOc, Gp and 9c being Fermi's cou- 
pling constant and Cabibbo's angle, Ei is the energy of 
the final state lepton and k and k' are the neutrino and 
charged lepton momenta, respectively. Compared to the 
corresponding quantities appearing in Eq. (0), the ten- 
sors L^^ and W^'^ include additional terms resulting from 
the presence of axial-vector components in the leptonic 
and hadronic currents (see, e.g., Ref. (59l]). 

Within the lA scheme, the cross section of Eq. H37|) 
can be cast in a form similar to that obtained for the case 
of electron- nucleus scattering (see Eq. ((T7|l ). Hence, its 
calculation requires the nuclear spectral function and the 
tensor describing the weak charged current interaction of 
a free nucleon, w'^ . In the case of quasi-elastic scat- 
tering, neglecting the contribution associated with the 
pseudoscalar form factor Fp, the latter can be written in 
terms of the nucleon Dirac and Pauli form factors Fi and 
F2, related to the measured electric and magnetic form 
factors Ge and Gm through 
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Ge) 



(38) 



(39) 



and the axial form factor Fa- 

FigurelTSlshows the calculated cross section of the pro- 
cess ^^0(ye, e), corresponding to neutrino energy E^ — 1 
GeV and electron scattering angle 9^ — 30° , plotted as a 
function of the energy transfer v = E^ — Ef,. Numerical 
results have been obtained using the spectral function 
of Fig. 121 and the dipole parametrization for the form 
factors, with an axial mass of 1.03 GeV. 

Comparison between the solid and dashed lines shows 
that the inclusion of FSI results in a sizable redistribu- 
tion of the lA strength, leading to a quenching of the 
quasi-elastic peak and to the enhancement of the tails. 
For reference, we also show the cross section predicted 
by the FG model with Fermi momentum pp = 225 MeV 
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FIG. 13: (Color online) Differential cross section da/dTledv 
for neutrino energy Ev = 1 GeV and electron scattering angle 
6e — 30° . The lA results are represented by the dashed line, 
while the solid line corresponds to the full calculation, includ- 
ing the effects of FSI. The dotted line shows the prediction 
of the FG model with Fermi momentum kp = 225 MeV and 
average separation energy e = 25 MeV. 



and average separation energy e = 25 MeV. Nuclear dy- 
namics, neglected in the oversimplified picture in terms 
of noninteracting nucleons, clearly appears to play a rel- 
evant role. 
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FIG. 14: (Color online) Differential cross section da/dQ^ for 
neutrino energy E = 1 GeV. The dot-dash line shows the I A 
results, while the solid and dashed lines have been obtained 
using the modified spectral function ef Eq. H4()|l . with and 
without inclusion of FSI, respectively. . 



It has to be pointed out that the approach described 
in Section ^1 while including dynamical correlations in 
the final state, does not take into account statistical cor- 
relations, leading to Pauli blocking of the phase space 
available to the knocked-out nucleon. 

A rather crude prescription to estimate the effect of 
Pauli blocking amounts to modifying the spectral func- 



tion through the replacement 

P{ip,E)^P{p,E)e{\p + q\-pp) (40) 

where pp is the average nuclear Fermi momentum, de- 
fined as 



Pf 



(fr pA{r)pF{r), 



(41) 



with p_F(r) = (37r^pA(r)/2)^/'^, being the nuclear 

density distribution. For oxygen, Eq. (|41|l yields pp = 
209 MeV. Note that, unlike the spectral function, the 
quantity defined in Eq. H40|l does not describe intrinsic 
properties of the target only, as it depends explicitely on 
the momentum transfer. 
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FIG. 15: (Color online) Quasi-elastic differential cross sec- 
tion da/dEf^ as a function of the scattered energy E^ for the 
neutrino energy E = 0.8 and 2.0 GeV. The solid line shows lA 
calculation with Pauli blocking as in Eq. (40), the dot-dash 
line lA calculation without Pauli blocking, and the dashed 
line FG model. 

The effect of Pauli blocking is hardly visible in the dif- 
ferential cross section shown in Fig. 1131 as the kinemat- 
ical setup corresponds to > 0.2 GeV^ at the quasi- 
elastic peak. The same is true for the electron scattering 
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cross sections discussed in the previous Section. On the 
other hand, this effect becomes very large at lower . 

Figure^lshows the calculated differential cross section 
da/dQ^ for neutrino energy — \ GeV. The dashed and 
dot-dash lines correspond to the lA results with and with- 
out inclusion of Pauli blocking, respectively. It clearly 
appears that the effect of Fermi statistic in suppressing 
scattering shows up at < 0.2 GeV^ and becomes very 
large at lower Q^. The results of the full calculation, 
in which dynamical FSI are also included, are displayed 
as a full line. The results of Fig. E| suggest that Pauli 
blocking and FSI may explain the deficit of the measured 
cross section at low with respect to the predictions of 
Monte Carlo simulations [Tll |. 

FigureElshows the i^^-nucleus cross sections as a func- 
tion of the scattered muon energy, by comparing the cross 
sections calculated by FG, and by the use of the spec- 
tral function with and without Pauli blocking. Figure ITKl 
shows that FG yields a larger high-energy peak contri- 
bution than the other two. This is not due to the PauH 
blocking, but due to the nuclear correlation effects in the 
spectral function: the muons tend to be scattered with 
a higher energy. This effect should show up in the for- 
ward angle cross section and may have a direct effect on 
neutrino oscillation measurements. 



V. CONCLUSIONS 

We have employed an approach based on NMBT to 
compute the inclusive electron- and neutrino-nucleus 
scattering cross sections in the kinematical region cor- 
responding to beam energy 1 GeV, relevant to many 
neutrino oscillation experiments. Our calculations have 
been carried out within the lA scheme, using realistic 
spectral functions obtained from (e, e'p) data and theo- 
retical calculations of uniform nuclear matter. 

In the region of the quasi-elastic peak, the results of our 
calculations account for the measured ^^0(e, e') cross sec- 
tions at beam energies between 700 MeV and 1200 MeV 
and scattering angle 32° with an accuracy better than 
10 %. It must be emphasized that the ability to yield 
quantitative predictions over a wide range of beam ener- 
gies is critical to the analysis of neutrino experiments, in 
which the energy of the incident neutrino is not known, 
and must be reconstructed from the kinematics of the 
outgoing lepton. 

In the region of quasi-free A production theoretical 
predictions significantly underestimate the data. Assum- 
ing the validity of the lA scheme, this problem appears 
to be mainly ascribable to uncertainties in the descrip- 
tion of the nucleon structure functions in this kinematical 



regime. The upcoming electron-nucleus scattetring data 
in the resonance region from the Jefferson Lab E04-001 
experiment will help to shed light on this issue. At 
higher energies, i.e. in the region in which inelastic con- 
tributions largely dominate, the calculated cross sections 
are in close agreement with the data. 

Among the mechanisms not included in the lA picture, 
scattering processes in which the incoming lepton cou- 
ples to meson exchange currents are not expected to pro- 
duce large corrections to our results in the region of the 
quasi-elastic peak. Numerical studies of the transverse 
response of uniform nuclear matter, carried out within 
NMBT jsJl , have shown that inclusion of two-body con- 
tributions to the nuclear electromagnetic current, arising 
from TT and p meson exchange, leads to an enhancement 
that decreases as the momentum transfer increases, and 
never exceeds 10 % at < o.25 GeV^. On the other 
hand, the results of calculations of the tranverse response 
of the few- nucleon systems (for a review see Ref. j2d| ) 
suggest that two-body current contributions may play a 
role in the dip region, at least for the lower values of the 
momentum transfer. 

The second mechanism not included in the lA, Pauli 
blocking, while not appreciably affecting the lepton en- 
ergy loss spectra, produces a large effect on the distri- 
butions at < 0.2 GeV^, and must therefore be taken 
into account. 

In conclusion, NMBT provides a fully consistent and 
computationally viable scheme to calculate the elec- 
troweak nuclear response. Using the approach discussed 
in this paper may greatly contribute to decrease the sys- 
tematic uncertainties associated with the analysis of neu- 
trino oscillation experiments, as, unlike the FG model 
and other many-body approaches based on effective NN 
interactions (see, e.g., Ref. [g^ and references therein), 
it is strongly constrained by NN data and involves no 
adjustable parameters. 
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